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Abstract. We investigate consequences of a continuously energy- injecting central engine of gamma-ray burst 
(GRB) afterglow emission, assuming that a highly magnetized pulsar is left beaming in the core of a GRB pro- 
genitor. Beaming and continuous energy-injection are natural consequences of the pulsar origin of GRB afterglows. 
Whereas previous studies have considered continuous energy-injection from a new-born pulsar to interpret the 
deviation of afterglow light curves of GRBs from those with the simple power law behavior, a beaming effect, 
which is one of the most important aspects of pulsar emissions, is ignored in earlier investigations. We explicitly 
include the beaming effect and consider a change of the beaming with time due to a dynamical evolution of a 
new-born pulsar. We show that the magnitude of the afterglow from this fireball indeed first decreases with time, 
subsequently rises, and declines again. One of the most peculiar optical afterglows light curve of GRB 970508 can 
be accounted for by continuous energy injection with beaming due to a highly magnetized new-born pulsar. We 
discuss implications on such observational evidence for a pulsar. 

Key words, gamma rays:bursts - pulsar:general ~ stars:magnetic fields 



1. Introduction 

. Gamma-ray bursts (GRBs) are widely accepted to be pro- 
I duced when fast-moving, relativistic shells ejected from 
1 a central source in a relatively short period collide with 
slowly moving, yet relativistic shells that were ejected at 
an earlier time (Rees & Meszaros 1994; Paczynski & Xu 
1994; Kobayashi et al. 1997; Daigne & Mochkovitch 1998, 
2000). In connection with the so-called internal shock 
model, the external shock model also prevails as a possible 
origin of the GRB afterglows. In the external shock model 
the relativistic material is assumed to be decelerated via 
interactions with the surrounding medium. A shock wave 
results in heating the ambient matter to relativistic tem- 
peratures, and emitting photons in longer wave lengths 
(Rees & Meszaros 1992; Meszaros & Rees 1993; Paczynski 
& Rhoads 1993; Sari et al. 1996; Vietri 1997). GRBs and 
their afterglows seem to result from the dissipation of bulk 
energy in the relativistic outflows, which are in the form 
of a narrow beam rather than a spherical shell. 

Even though the origin of the observed GRBs are still 
unknown, from the observations of several GRB after- 
glows the evidence of beamed GRBs has been accumu- 
lated (Sari et al. 1999; Halpern et al. 1999; Rhoads 1999). 
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There are several works on models for the geometry of 
GRBs (e.g., Chang & Yi 2001 and references therein) and 
their environments (e.g., Scalo & Wheeler 2001 and ref- 
erences therein). Much of the current research on GRBs 
is aimed at determining the nature and the origin of the 
central engine (Duncan & Thompson 1992; Narayan et al. 
1992; Woosley 1993; Katz 1994; Usov 1992, 1994a; Shaviv 
& Dar 1995; Meszaros & Rees 1997b; Yi & Blackman 
1997; Blackman & Yi 1998; Paczynski 1998; MacFadyen 
& Woosley 1999; Portegies Zwart et al. 1999; Li 2000; 
Wheeler et al. 2000; Zhang & Fryer 2001). Although the 
simple cosmological fireball afterglow model is in a good 
agreement with the observed light curves of a power law 
decay (e.g., Wijers et al. 1997; Waxman 1997a, b), the af- 
terglow of GRB970228 observed by HST, for instance, 
shows a deviation from the simple power law behavior 
(Fruchter et al. 1997). Several works have been done to 
further investigate more subtle effects that can change the 
afterglow characteristics (Katz & Piran 1997; Meszaros & 
Rees 1997a, 1999; Rhoads 1999; Berger et al. 2000; Dai 
& Lu 2000; Kumar & Panaitescu 2000a, b; Panaitescu & 
Kumar 2000; Dai & Lu 2001). 

Most fireball models assume that the energy injection 
into the fireball occurs in a short period of time compared 
with the lifetime of the afterglows (Rees & Meszaros 1998; 
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Kumar & Piran 2000; Sari & Meszaros 2000). However, in 
some types of central engines, such as a fast rotating new- 
born pulsar with the strong magnetic field (magnetar), a 
significant energy input into the fireball may in principle 
continue for a significantly longer timescale, and accord- 
ingly the temporal decay of the afterglow will be slower. 
Hence, it is worthwhile to investigate a continuously fed 
fireball in more details as a probe of the central engine 
of GRBs. We here consider the central engine that emits 
both an initial impulsive energy input i^imp and a con- 
tinuous power. In fact, recently there were such attempts 
to provide an explanation for the deviation of the after- 
glow light curve from the simple power law (Dai & Lu 
1998a, b; Zhang & Meszaros 2001a, b). Even though pre- 
vious studios have considered continuous injection from a 
highly magnetized millisecond pulsar to interpret the de- 
viation of afterglow light curves of some GRBs from the 
power law behavior, one of the most important aspects 
of pulsar emissions is ignored, that is, beaming. Beaming 
and continuous powering are the clear consequences of the 
pulsar origin of the GRB afterglows. In this Letter, this 
feature is included and a change of the beaming due to a 
dynamical evolution of the new born pulsar is explicitly 
taken into account. We find that the deviations suggested 
by previous studies indeed occur, but that the shape of 
the light curve is significantly modified due to a beam- 
ing effect. If such observational evidence for beaming is 
found, the corresponding pulsar origin for GRBs would 
imply that activities resulting in too massive objects to 
be a neutron star might be ruled out as a central engine 
of GRBs. Therefore, the question of whether or not a sig- 
nature of beaming can be found in the afterglows thus 
extremely important. 

2. Afterglows Continuously Powered by Pulsar 

The total luminosity emitted from a young millisecond 
pulsar (MSP) has two important terms: an electromag- 
netic (EM) radiation term and a gravitational wave (GW) 
radiation term. Given that the spin-down is mainly due 
to electromagnetic dipolar radiation and to gravitational 
wave radiation, the spin-down law is given by 
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where H, and Cl are the angular frequency and its time 
derivative, Bp is the dipolar field strength at the mag- 
netic poles, a is the angle between spin axis and the 
magnetic field axis, / is the moment of inertia, R is 
the pulsar radius, and e is the ellipticity of the com- 
pact star (Shapiro & Teukolsky 1983). When an MSP 
rotates with spin frequency larger than a certain criti- 
cal frequency, flait, nonaxisymmetric secular instability 
drives the pulsar into nonaxisymmetric configuration with 
nonzero quadrupole moment. This is the case where the 
GW luminosity dominates. In this case, the spin-down 



of the MSP occurs on a timescale r. 
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MSP moment of inertia, = f24lO^ s ^ is the MSP spin 
frequency, and we have made use of the gravitational wave 
energy-loss rate, Lgw = 32Ge'^linl/5c^ . When an MSP 
rotates with spin frequency smaller than the certain crit- 
ical frequency, the spin-down due to the electromagnetic 
dipole radiation gives a timescale Tem ~ I^fl1/2Lcm ~ 
3 X lO'^h^B^^^R^^n^'^ s, where i?* = RqIO^ cm is the 
MSP radius, B^ = Bi^lO^^ G is the MSP dipole magnetic 
field, and we have assumed the electromagnetic dipole 
energy-loss rate, Lem = 2/i2ri^/3c'^, where /i* = B^R^ 
is the electromagnetic dipole moment. The decay solution 
0(T) as a function of time T includes both EM and GW 
emissions, but the energy input into the fireball is due 
to the EM dipolar emission Lem alone. Provided that the 
spin-down at a time is dominated by one or the other 
loss term, one can obtain approximate solutions. When 
GW radiation losses dominate the spin-down, we have 
Q = f2o,gw(l + T/Tgw)~^/^, or approximately fl = Oo.gw 
for T < Tgw, and ^ = ^^o.gwCT/rgw)"^/" for T » Tgw, 
where f^o.gw being the initial spin frequency. When EM 
dipolar radiation losses dominate the spin-down, we have 

Q = f2o,em(l + T/Tem)^^^^ , OT approximately fl = flo,eni 
for T < Tem, and O = r2o,em(T/Tem)"^/^ for T » Tem, 

where ^o,em being the initial spin frequency. Assuming 
i^crit > ^^dynamo, where ^dynamo is the Spin required for 
dynamo action, there are two types of MSPs: a supercrit- 
ical strong field rotator (SPS) with O > ftcHt > i^dynamo, 
and a subcritical strong field rotator (SBS) with ficrit > 
SI > ridynamo- Thcsc two classcs may give the bimodal 
distribution of short and long pulsar-induced GRBs (Yi 
& Blackman 1998). For an SPS pulsar, GW spin-down is 
important in initial period and promptly below flcrit- As 
soon as < flait the spin-down becomes dominated by 
the EM losses. On the other hand, for an SBS pulsar, the 
spin-down is always dominated by EM regime. 

The differential energy conservation relation for the 
self-similar blast wave can be written as dE/dt = 
Coit/taY — i^'{E/t), where E and t are the energy and 
time measured in the fixed frame and q' and k' are con- 
stants (Cohen & Piran 1999). The first term denotes 
the continuous luminosity injection, and the second term 
takes into account radiative energy losses in the blast 
wave. For t > Iq, the bulk Lorentz factor of the fire- 
ball scales with time as oc with m and k,' re- 
lated by k' = m — 3 (Cohen et al. 1998). If m = 3, it 
corresponds to the adiabatic case (Blandford & McKee 
1976). In the observer frame, the time T is related to 
the fixed frame t hy dT = (1 - P)dt ~ dt/2V'^, and 
T = J^{2T^)-^dt ~ t/[2{m + l)r^] when t > to- The 
differential energy conservation relation in the observer 
frame is now given by dE/dT = LoiT/To)" - k{E/T), 
and can be integrated as 



E 



tY [t 

— \T+ — 



(2) 



gw 



3 X 10-3e-2/45f]~'^ s, where J* = /45l0''^ g cm^ is the 



where Tq — io/[2(?Ti' + l)r^], m' being the self-similar 
index for t < to, Lq = 2r^£o, Q = {q' — m)/{m + 1), and 
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K = K,' /{m + 1). Setting T = Tq, the total energy at the 
beginning of the self-similar expansion is the sum of two 
terms, Eq = LqTq / (k + q + 1) + Ei^np. The first term is the 
accumulated energy from the continuous injection before 
the self-similar solution begins. The second term i?imp is 
the energy injected impulsively by the initial blast. 

The total energy of the blast wave given by Eq. ^ 
may be dominated either by the continuous injection term 
(oc yj+i) or by the initial impulsive term (oc T^^), sub- 
ject both to the relative values of the two indices and to 
the values of Lq and i?imp (Dai & Lu 1998a, b; Zhang & 
Meszaros 2001a). One may classify three regimes accord- 
ing to the relative values of the two indices as discussed by 
Zhang & Meszaros (2001a). We are interested in the case 
where q > — 1 — k, since otherwise a pulsar signature is no 
longer observable even if there is a pulsar in the central 
engine. If g > — 1 — k, the first term in Eq. |^ will even- 
tually dominate over the second term after a critical Tc- 
The injection-dominated regime begins at a critical time 
Tc defined by equating the injection and energy-loss terms 
in Eq. I 



Tr = Max<^ 1 



(K + q + l) 



Es 



imp 



l/{K+q+l) 
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If initially the impulsive term dominates, the critical time 
Tc could be much longer than Tq, depending on the ra- 
tio of -Eimp and LqTq. The case where > Tq ensures 
that a self-similar solution has already formed when the 
continuous injection term dominates. Besides the condi- 
tion oi q > — 1 — K, for detecting the signature of a pulsar 
the critical time Tc should be greater than Tom or Tg^ as 
discussed below. 

To obtain the temporal decay index of the afterglow 
light curve for which the MSP is responsible we adopt 
the cylindrical geometry instead of the spherical geome- 
try, which may accommodate elongated beaming config- 
urations. The rotational axis of the MSP coincides with 
the z-axis of the geometry. For a fireball blastwave decel- 
erated by a homogeneous external medium with particle 
number density n, the energy conservation equation at 
time t = r/c is given by 



E = 2tt R^rnmpC^V'^ = 271 ctnmpC^T'^ , t > to, 



(4) 



where r is the distance of the blastwave in z direction, 
i? is a radius of the beaming which can be characterized 
by that of the light cylinder, c/il, being given by the 
MSP evolution, and all other symbols have their usual 
meanings. One can work out the relationship between 
the temporal index a and the spectral index /3, where 
F„ cx T"!^^. To directly compare indices suggested by 
Zhang & Meszaros (2001a), we derived the scaling laws for 
the slow cooling regime assuming that the reverse shock 
is mildly relativistic (Sari et al. 1998). We derive the scal- 
ing law for T Tom as an illustrating example. The en- 
ergy E in Eq. ^ should have the same time dependence 
as the first term on the right-hand side of Eq. |^, giv- 
ing T'^+^ (X iP^, or TO = —q/{q -\- 2). Since in general, 



Table 1. Temporal index a for the forward and the reverse 
shocks in various timescales. 











T < Tgw 


3-f2q 




(9g + 12 - 2g/3'-)/4 


T » Tgw 


2 + 2g 




(18g -H5 - 2(2g - l)/3'')/8 


T < Tcm 


3 + 2g 


~ql3t 


(9g + 12- 2g/3' )/4 


T > Tcm 


l + 2q 


-pf{q^l) 


(9g + 3-2/3'-(g-l))/4 



P dependence on T is given by P (X T'/^. The relationship 
between the temporal index a and the spectral index jS 
can be obtained as below. For the forward shock, the syn- 
chrotron peak frequency f/„ cx P* cx r-2m/(m+i) qc T^, 
the peak flux F^^ cx T^T^ cx T'^3-m)/{m+i) ^ j^3+2q^ 
Thus, af ^ 3 + 2q - q(3 = {2mf3-f + 3 - m)/(l -f to) 



since T' 



y3+2g ^ y(3-m)/(m+l)^ p^^, rCVerSC 



shock, ly:;^ = vlJT^ cx P^ cx y-^/^^+i) cx T^/^, the peak 

flux = PF/^ cx r3p9 (X r(6-3™)/2(m+l) T(l2+9q)/i^ 

Thus, a" = (2TO'i9''-h6-3TO)/2(l-t-TO) = (12 + 9g-2(j/3)/4. 
Similarly, one can derive the temporal index for each 
regime. In Table 1 we have summarized temporal indices 
in various regimes for the forward shock and the reverse 
shock, respectively. 

If the continuous injection term becomes dominant 
over the impulsive term after Tc, the afterglow light curves 
rises after Tc and steepen again after some time, that is, 
about Tom or Tg^. Therefore, there may be two types of 
afterglow patterns for continuous injections according to a 
appropriate combination. Zhang & Meszaros (2001a) dis- 
cussed conditions which allow to detect a signature for a 
pulsar and concluded that physical parameters are con- 
sistent with those of a magnetar in case of the afterglow 
features in case of GRB 000301c. In practice, however, Tgw 
is very short, can be even shorter than Tc unless ambient 
matter density is very high, and is therefore unlikely to be 
observed. At around Tc and Tom, ? changes -1 to and 
to -2, respectively. These scaling laws represent a change 
from the standard adiabatic case to an EM-loss dominated 
regime as shown above (see also Zhang & Meszaros 2001a). 
For the forward shock, the temporal decay index changes 
around Tc from 3(3^/2 to 3 and returns to 3f3-^ — 3 af- 
ter Tom. For the reverse shock, the temporal decay index 
changes around Tc from 3/3''/4 — 3/8 to 3 and returns to 
3/372 - 15/4 after Tom. 

3. Discussions 

A continuous energy injection signature in the GRB af- 
terglow light curve may directly provide diagnostics about 
the nature of the injection as well as information on the 
GRB progenitor. Therefore, the question of whether or 
not a bump in the afterglow light curve is such an ob- 
servational evidence for beaming from a pulsar is thus 
extremely important. In this sense, the correct geometry 
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should be applied to the model calculation. We have dis- 
cussed the case where a pulsar is continuously injecting 
energy cylindrically rather than spherically. We suggest 
that a possible explanation for the deviation of the after- 
glow light curve of GRB 970508 is due to a beaming from 
a central pulsar. The optical afterglow of GRB 970508 has 
been explained as evidence of gravitational lensing event 
(Loeb & Perna 1998; Dado et al. 2001). So far the gravita- 
tional lensing can only account for the fact that the GRB 
afterglow shows a rise and a fall. We show in this Letter 
that with a correctly assumed geometry of energy injec- 
tion the afterglow of GRB 970508 can be explained by 
a strongly magnetized fast-rotating pulsar with the con- 
tinuous energy injection in a beam. In this calculation 
we ignore effects of the lateral expansion of the jet which 
may occur in the later time of the jet evolution. The light 
curve should be modified and becomes that of the spher- 
ical geometry case if the sideways expansion occurs in a 
timescale comparable to those we discussed after a neu- 
tron star formed. For instance, if the jet is expanded much 
faster at earlier stage than the light cone evolution, the 
cylindrical geometry effect becomes less obvious. 

If the bump is caused by such a pulsar indeed, it puts 
constraints on the GRB progenitor models. That is, the 
corresponding pulsar origin for GRBs would imply that 
activities resulting in too massive objects to be a neu- 
tron star, such as, neutron star mergers, black hole for- 
mation, should be ruled out as a central engine of GRBs. 
During the birth of the neutron star, an initial fireball 
may occur through electromagnetic processes (Usov 1992). 
This has led to models in which GRBs are powered by 
rapidly spinning compact objects with strong magnetic 
fields (Blackman & Yi 1998; Blackman et al. 1996; Usov 
1994b; Yi & Blackman 1997, 1998). 
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